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Abstract. In this paper we study the Cauchy problem for an p— Laplacian type 
of evolution system H 4 + V x [|V x H| P ~ 2 V x H|] = F. This system governs 
the evolution of a magnetic field H, where the current displacement is neglected 
and the electrical resistivity is assumed to be some power of the current density. 
The existence, uniqueness and regularity of solutions to the system are established. 
Furthermore, it is shown that the limit solution as the power p — > oo solves the 
problem of Bean's model in the type-II superconductivity theory. The result pro- 
vides us information about how the superconductor material under the external 
force to become the normal conductor and vice visa. It also provides an effective 
method to find numerical solutions to Bean's model. 



1. Introduction 
In this paper we shall study the following degenerate evolution system: 

H t + V x [| V x H| P ~ 2 V x H] = F, (x, t) G Q T , (1.1) 
V-H = 0, (x,t)eQ T , (1.2) 

H(x,0) = H (ar), x G R 3 , (1.3) 

where Qt = R 3 x (0, T] for some T > and p > 2. Here a bold letter represents a 
vector in R 3 . 

The system (1.1) is derived from Maxwell's system where the current displacement 
is neglected since it is small in comparison of eddy currents ( see Landau- Lifschitz 
[12]). The vector H represents the magnetic field. The electrical resistivity is assumed 
to be equal to some power of the current density. The motivation for our investigation 
of (1.1)-(1.3) is twofold. On one hand, the system (1.1) is a natural generalization 
of the scalar p-Laplacian which has been studied by many authors (see DiBenedetto 
[6] and the references therein). More recently, the limit solution of the p-Laplacian 
equation as p — > oo is used to describe the fast-slow diffusion model and the sand- 
pile collapsing phenomenon (see Arronson-Evans-Wu [1], Evans- Gangbo [9], Evans- 
Feldman-Gariepy [8] and the references therein). On the other hand, it will be seen 
that the limit solution to the system (1.1)-(1.3) solves the problem of Bean's critical- 
state model for type-II superconductor (see Bean [3]). In Bean's model, the current 
density is assumed to be always less than or equal to a critical current, denoted by 
J c . The superconductive region is characterized by the current density J = |V x H|. 
When the current density is strictly less than the critical number J c , then there is 
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no resist for the movement of electrons (the superconductivity is achieved) while the 
normal conductor region is the one where \J\ = J c . When a magnetic field is under 
the influence of an external force, a normal conductor material will become a super- 
conductor one or vice visa. The motion of the interface is driven by the external force. 
Recently, Prigozhin in [14] reformulated the problem to a variational inequality and 
established the existence of a unique weak solution to the variational inequality if 
the space dimension is equal to two. Some numerical issues are discussed in [2] (also 
see a more recent paper [15]). Surprisingly, we show that the weak solution defined 
in [14] can be obtained as the limit of the solution to the system (1.1) as p — > +oo. 
Here we assume that the Bean's critical number J c = 1 without loss of generality. 
Our argument is entirely different from that of [14]. From the physical point of view, 
for large p the electric resistivity is small in the region S e — {(x, t) : | V x H| < 1 — e} 
while it is very large in {(x, i) : | V x H| > 1 + e}. Thus, the resistivity in S e becomes 
smaller and smaller as p increases. And S £ becomes the superconductor region as 
e — > (no resistivity). The region {(x,t) : 1 — e < |V x H| < 1 + e} plays as the 
intermediate zone and the limit of the zone as e — > becomes the interface bewteen 
the normal and superconductor regions. 

Unlike the scalar p-Laplacian equation there is not much work being done about the 
degenerate system (1.1) as well as its limit problem. Recently, the regularity of weak 
solutions to a linear system analogy to (1.1) was studied in [17]. In this paper we shall 
first prove the existence of a weak solution to the system (1. l)-( 1.3) . Then we study 
the limit solution as p — > oo and prove that the limit solution solves a variational 
inequality ([14]). Many techniques are adopted from the scalar p-Laplacian case ([8]). 
To study the regularity and understand the profile of the limit solution, we consider 
the magnetic field H to be plane waves. Then it is rather surprising to see that the 
current density satisfies a porous medium equation. The limit problem for (1.1)- (1.3) 
becomes the mesa problem ( See Caffarelli- Friedman [5]) for the density equation. 
This allows us to analyze the structure of the solution to the limit problem. 

The paper is organized as follows. In § 2, by using standard variational method 
we show that there exists a unique weak solution to the system (1.1) -(1-3). Then 
we prove that the current density satisfies a porous-medium type equation when the 
magnetic field is assumed to be plane waves. For this special case, further regularity 
of the weak solution is obtained. In § 3, we prove that the solution of (1.1)-(1.3) has 
a unique limit as p goes to infinity. Moreover, the limit solution solves a variational 
inequality. In § 4, we study the mesa problem with an inhomogeneous term. Similar 
results to the paper [5] are obtained. In § 5, we use the results obtained in §4 to 
study the structure of the limit solution of (1.1)-(1.3). 

2. Existence, Uniqueness and Regularity for fixed p > 2 

Let p > 2. Define 

B d {R 3 ) = {g(x) e W^iR 3 ) : V • G = 0, a.e.x e R 3 } . 

We shall assume the following conditions on H and F. 
H(2.1) Assume that H e B d (R 3 ),F e L 2 {0,T; B d (R 3 )). 
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Definition 2.1: A vector field H G L 2 (0,T; B d (R 3 )) is said to be a weak solution of 
the problem (l.l)-(l.S), if the following integral identity holds: 

£ J [-H • $t + | V x H| p - 2 (V x H) • (V x $)] dxdt 

= J [H (a;) • $] rfx (2.1) 

/or a// *) G H\0, T; B d (R 3 )) with $(x, T) = on R 3 . 
First of all, we derive some elementary energy estimates. 

Lemma 2.1: Under the assumption H(2.1) there exist constants C\ and C 2 such 
that 



sup 

ro,n J r 3 



£ \H\ 2 dx + J J q |V x H\ p dxdt < Ci, (2.2) 

// |H| 2 (M+ isup / IV x H| p (ix < C 2 . (2.3) 
J JQt V fn.n JR' A 



V [0,T] 

where C\ and C2 depend only on known data and the upper bound of T . 
Proof: Note that for any vector fields A, B G B d the following identity holds: 



|^A-(Vx B)dx = ^B-(Vx A)dx. 



IR 3 J R? 

Taking inner product to the system (1.1) by H and then using the above identity, we 
obtain 



sup / |H| 2 cfe + [ T [ |V x B.\ p dxdt 
[o,t] JR 3 Jo Jr 3 

< [ \H \ 2 dx+ [ T [ F-Udxdt. 

Jr 3 Jo Jr 3 

By first using Cauchy's inequality and then using Gronwall's inequality, we have 

sup / \H\ 2 dx + [ T ( IV x B.\ p dxdt 
[o,t] Jr 3 Jo Jr 3 



< [ |H | 2 dx + C ( T [ \F\ 2 dxdt 
Jr 3 Jo Jr 3 



To obtain the second estimate, we take the inner product by H t to (1.1). Similar to 
the first estimate one can easily derive the second estimate. 

Q.E. D. 

Theorem 2.2: Under the assumption H(2.1) the problem (l.l)-(l.S) has a unique 
weak solution. Moreover, 

H t G L 2 (Q T ),H G L°°(0, T; B d (R 3 )). 
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Proof: The proof is based on the standard variational method. For every V G 
L 2 (0,T; R 3 ), define 



Ip\Y\ 



lf R3 \VxV\Pdx,UVeB d , 
+00, otherwise, 



for a.e. t G (0,T]. The system (1.1) can be reformulated by the following variational 
inequality: 

F - Hi G dl p [H] , x G Z? 3 ,fora.ei G [0,T], (2.4) 
H(x,0) = H , xGi? 3 , (2.5) 

where <9/ p [H] is the sub differential of J P [V] at V (see [16]). 

It is clear that the operator I p [V] is coercive and strictly monotonic increasing. 
Moreover, The operator I p is mapping from L 2 (0,T; W 1 ' P (R 3 )) to 

L 2 (0,T; W^'P'iR 3 )), where p' = It follows ([16]) that the evolution problem 
has a unique weak solution H G L 2 (0,T; B d ). Moreover, by Lemma 2.1 the weak 
solution H t G L 2 (Q T ) and V x H G L p (R 3 ) for a.e. t G [0, T\. 

Q.E.D. 

Further regularity of the weak solution seems challenging since there are no local 
estimates such as those in [6]. However, we can obtain more regularity when the 
weak solution H is assumed to be plane waves. 

From now on we shall assume that H depends only on x = (xi,x 2 ) and the 
component in z-direction is zero, i.e., H(x,t) = {h±(x, t), fi2(x, t), 0}. In this case 
Qt = R 2 x{0,T]. 

To obtain more regularity of the weak solution, we need more regularity for the known 
data. 

H(2.2): Assume that H G C 2+a (R 2 ) and F G C 2+a (Q T ) for some a G (0, 1). 

Theorem 2.3: Under the assumptions H(2.1)-H(2.2), the weak solution of the prob- 
lem (l.l)-(l.S) belongs to L°°(Q T ) f| C 1+Q?,/3 (Q T ) for some a,p G (0,1). Moreover, 
the solution has compact support for each t G [0,T] ; ifF and H have compact sup- 
port. 

Proof: Recall that H = {hi(x, £), h 2 (x, t), 0}. Then it is clear that 

V x H = (h 2xi - h 1X2 )k, 
where k is the unit vector in z-direction. It follows that 
V x (|V x H| P ~ 2 V x H) 

— \ \h-2x! — hi X2 \ p 2 (h 2xi — h lx2 ) ,— \h 2xi — h lx2 \ p 2 {h 2xi — h lx2 ) ,0> 
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ox 2 

hot - 



\h 2xi 


-h 1X2 r 


2 (h 2xi 


— hi X2 ) 


= fi(xi 


x 2 ,t), 


(2.7) 


\h 2x x 


-h 1X2 r 


2 (h2 Xl 


— hi X2 ) 


= hixx 


x 2 ,t). 


(2.8) 



From Theorem 2.2, we know that hi(—oo, x 2 , t) — for all x 2 G R l and t > 0. Then 
from the equation (1.2), we see 

/xx 
h 2x2 (i,x 2 )di = -U 2x2 (x 1 ,x 2 ,t), 
-co 

where 

/xx 
h 2 (£,x 2 ,t)d£. 
-co 

It follows that 

hi X2 = —U 2x2X2 ,andh 2xi = U 2xixi ,(xi,x 2 ,t) G Qt- 
Consequently, the equation (2.7) becomes 

d 



dx 2 



-U 2t + \AU 2 \ p - 2 AU 2 =f 1 (x 1 ,x 2 ,t). 



Integrating over (— oo,x 2 ) yields 

U 2t -\AU 2 \ p - 2 AU2 = f;(x,t), 
where 

&{x,t) = - H f 1 (x 1 ,U)dC 

J —CO 

Similarly, we can eliminate h 2 (x,t) to derive the following equation from (2.8): 



(2.9) 



U lt -\MJ 1 \*~ 2 MJ 1 = ft(x,y,t), 



(2.10) 



where 



/XI 
hi(xi,Z,t)d£, 
-co 

K(x,t) = - [ X1 m,x 2 ,t)dt. 

J —CO 

Define 

m(x,t) = AUi{x,t), 9i (x,t) = Af*(x,t), i = l,2. 
Then U\ and u 2 will satisfy the following system: 

u it - A^m) = gi(x,t), (x, t) G Qt, 

Ui(x, y,0) = u i(x), xeR 2 , i = l,2, 

where 

/xi rxx 
h 01 (x l7 ^)d^,u 02 (x) = A h 02 (C,x 2 )d^ 
-CO J —CO 

/x\ rX2 
f2^,X2,m,g 2 (x,t) = -A / fi(x u €,t)d€, 
-CO J —CO 



(2.11) 
(2.12) 



9i ( 
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while the function tp is defined by 

V>(s) = |s| p_2 s. 

By a result of [7] (also see [10]), there exists a unique solution 

Ui(x,t) eL 2 (Q T )f)L°°(Q T ). 

Moreover, the solution is Holder continuous in Qt- Furthermore, the suppUi is com- 
pact for each t G [0, T], if ^ and -Uo« have compact support, i = 1,2. 
Now for each fixed t G (0, T], we solve the following elliptic problem: 

AUi = u h x G R 2 , (2.13) 
Ui(x,t)=0, \x\ -> +oo. (2.14) 

By applying the standard elliptic theory [11], we know that the problem (2.13)-(2.14) 
admits a unique solution Ui(x,t) G C 2+a '°(i? 2 ). Now we claim that Ui(x,t) is also 
Holder continuous with respect to t. Indeed, for small At we define 

Ui(x,t + At) -Ui(x,t) 



U* = 



Then U* solves the following problem: 



\At\ a / 2 



AU* = ^±^^ ,xeR 2 , (2.15) 
\At\ 2 

U*(x,t) = 0, |x|^+oo. (2.16) 

As Ui(x,t) is Holder continuous with respect to t, we apply L p -theory for elliptic 
equations to obtain: 

H^llwKB^o)) < C, 

where B R (0) = {x G R 2 : \x\ < R} and the constant C depends only on p, Holder 
norm of Ui and R, but independent of At. 

By Sobolev's embedding, we know that for every R > there exist constants C(R) 
and a G (0, 1) such that 

II^IIci+^b^o)) < C(R). 

A compactness argument yields that Ui(x,t) and Ui Xj are Holder continuous with 
respect to t for each i and j, where % — 1,2 and j = 1,2. It follows that hi(x,t) = 
-r^U\{x,t) and h 2 = -£^U2(x,t) are Holder continuous with respect to t. Finally, an 
interpolation lemma ([13]) implies that hi is Holder continuous with respect to t with 
Holder exponent (3 for some (3 G (0, 1). 

Q.E.D. 

3. The Limit Solution as p — > oo 

In this section we shall show that the solution of (1.1)-(1.3) has a limit as p — > oo, 
which solves a variational inequality if the initial current density is less than the Bean 
critical value. 
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H(3.1): Assume that the initial field H satisfies 

II V x H ||L°°(i?3) < 1. 

Lemma 3.1: Under the assumptions H(2.1) and H(3.1), the energy estimates (2.2)- 
(2.3) hold, where the constants C\ and C2 are independent of p. 

The proof is the same as for Lemma 2.1 except that the assumption H(3.1) will 
ensure that C\ and C2 are independent of p. 

We shall denote by the solution of the system (1.1)-(1.3). 

From Lemma 3.1, we know that every sequence H^ p •* with p' — > +00 has a subse- 
quence H( p ") which converges to a limit, denoted by H(°°\ weakly in * — L°°(f2x (0, T]) 
for any bounded domain f2 C i? 3 . Moreover, 

V x H (p,,) -> V x H (oo) , 

Hf > - Hj 00 ) 

weakly in L 2 (fi x (0,T)). Furthermore, since V • H = 0, H (p) e L 2 (0,T;B d ) and 
H t G L 2 (Q T ), it follows that 

H(p") -> H (oo) a.e. in Q T . 



Lemma 3.2: The limit H^ 00 ) satisfies the following estimate: 

ess sup I V x H (oo) | < I. 

Qt 

Proof: From the proof of Lemma 2.1, we see 

[ T [ |V x U\ p dxdt < Cp, 

JO JB? 

where C is a constant which is independent of p. 
Similar to [8], for every small 8 > we define 



Then, 



A s = {(x,t) € Qt : |V x H (p) | > 1 + 5} . 

(l + S)\A s \< J jf \VxH (p) \dxdt 

< lim inf / / IV x H_ Pt \dxdt 

< ^m^inf jf |V x H^l^dard^" l^] 1- ^ 

< |M 



It follows that 

|^|=0. 

The desired result follows since S is arbitrary. 



Q.E.D. 
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Before stating the main result in this section, we define another functional. Let 

K = {VGL 2 (Q T );V-V = 0}. 

For every V G K, define I P [V] as in Section 2 for a.e. t G (0,T]. Similarly, we define 

' 0, if V G K, |V x V| < 1, 
+oo, otherwise. 



/oo[V] 



Theorem 3.3: Under the assumptions H(2.1) and H(3.1) the limit function H^ 00 ) is 
unique and satisfies the following variational problem: 

F — H t G dloo [H] , fora.e.t>0, (3.1) 

H(x,0) = Ho(:c), xeR 3 , (3.2) 

that is, for a.e. t G [0,T], 

/ooM-JootH^)] (3.2) 
> J (F-H^) ■ (V-H^)dx a.e.t>0, (3.3) 

for all V G ^ with V t G L 2 (Q T ). 

Proof: Since the limit is unique we may assume the whole sequence H^ p ^ converges 
to H(°°). For simplicity, we denote by H the limit solution. Note that the system 
(1.1) is equivalent to the following variational form: 

F — H t G dip [H] . (3.4) 

By the definition of /^[V], we only need to verify the inequality for all V G 
L 2 (Q; R 2 ) with |V x V| < 1 a.e. t G [0,T]. On the other hand, Lemma 3.2 implies 
that Joo[H] = for a.e. t G [0, T\. Thus, the inequality (3.3) is equivalent to 

/ (F - H t ) • (V - H) dx < 0. 

J Ft? 

As solves the variational problem (2.4)-(2.5), we see 

I P [V] > / P [HW] + j R2 (F - H| p) ) • (V - H«) dx 

> J r3 (F - H^) • (V - H p ) dx. (3.5) 
By assumptions on V, we have 

lim JJV] = 0. 

On the other hand, since — ► H a.e and H G L 2 (Q T ), it follows that 

J r3 (F - H| p) ) • (V - H (p) ) dx -> J r3 (F - H t ) • (V - H) dx. 

After taking limit in the inequality (3.5), we obtain the desired inequality (3.3). To 
verify the initial condition, we note that 

/ \H(x,t) -n (x)\ 2 dx < cVt^o 
Jr 3 

as t -> 0, since ||H t || L2( Q T) < C. 
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Now we prove the uniqueness. Suppose there are two solutions H and H* to the 
limit problem. Then by Lemma 3.2 we know that 

|V x H| < 1, |V x H*| < 1. 

We choose H and H* as test functions in (3.3), respectively, to obtain 

J (F-H t ) ■ (y-H)dx < 0, 
F - H* t ) • (V - H.*)dx < 0, 



IR3 

Consequently, one obtains 



IR3 

It follows that 



/ (H* - H) t • (H - H*)dx < a.e. t > 0. 

J R 3 

' '|H - H*|| L 2 (Qt ) < 0, 



dV 

which implies the uniqueness. 

Q.E.D. 

Remark 3.1: The solution obtained in [14] is essentially the same as the limit 
solution in Theorem 3.3 if the space dimension is equal to 2. 

Now we investigate further regularity of the limit solution when H is a plane wave. 

Theorem 3.4: If the space dimension is equal to 2, then the limit solution H^ 00 ) is 
globally bounded and Holder continuous in Qt- 

Proof: When = {hi, h 2 ,0}, we can study the regularity for ■u-°°' ) . Recall that 
for each t G [0, T], 

AU^ P) =u'f\x,t), xeR\ (3.6) 
U?\x,t) =0, |x|->oo, (3.7) 

Since Hw^Hl^Qt) < 1,, the elliptic theory implies that for any q > 1 

\\Ui P \\w 2 'i{B R (0)) < C, 

where C depends only on R and q, but not on p. It follows by Sobolev's embedding 
and a compactness argument that for any a G (0, 1) U < f°\x,t) G C l+I3 '°(Q T ) for any 
(3 G (0, 1) and each t G [0,T]. It follows from the definition that G C Q '°(Q T ). 

Next we derive the regularity of U-°°\x,t) with respect to t. Recall that 
satisfies Eq.(3.6). Note that 

IKt||z,i(Q T ) < c, 
where C is independent of p. It follows that for any q G (1,2) 

\Wit ) \\w 1 'i(Q T ) ^ C> 

where C is independent of p. As U*- 00 -* G C 1+/3 '°(Qt), it follows by the interpolation 
result ([13]) that the limit function U-°°^ and U-^ are Holder continuous with respect 
to t. 

Q.E.D. 
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4. The Mesa Problem with an inhomogeneous Term 

As we have seen in section 2, the limit problem for the system (1.1)-(1.3) be- 
comes the mesa problem for the current density if the space dimension is equal to 
2. Throughout this section we shall study the mesa problem with an inhomogeneous 
term. By employing similar techniques to [5], we derive the profile of the limit solu- 
tion to the mesa problem. The result will be used to understand the initial collapsing 
and afterwards evolution process for the current density of the limit solution iT 00 ) in 
the case of two space dimensions (see Section 5 below). 

Consider the mesa problem with an inhomogeneous term. 

u t - Au m = g(x,t), (x, t) G R n x (0, oo), (4.1) 
u(x,0) = f(x), xeR n , (4.2) 

H(4.1). Assume that f(x) G L 1 (R n ) f] L°°(R n ) and 

g(x,t) G L\R n x (0,oo)) f]L°°(R n x (0,oo)). 
Moreover, f(x) > and g(x, t) > on Q T = R 2 x (0, T]. 

Lemma 4.1: Let 

M = ||/||l~(.r«) +T\\g\\ L ^ (QT) . 

Then 

\\u^\\ L oo {Qt) <M, 
Moreover, if M < 1, then has a unique limit u^°°\x,t) as m — > oo and 

u^\x,t) = f(x) + f g(x,T)dr, (x,t) e Q T . 
Jo 

Proof: As a first step, we may assume that f(x) and g(x, t) have compact support. 
Then u^ m \x,t) has compact support for each t G [0,T] and any m > 1. Let R be 
sufficiently large such that t) = on dB R (0) x [0,T]. We may also assume 

that 

u {m) is 

smooth. Define the operator L as follows: 

L[u\ = mu^Au + m(m - l)w m_2 |VM| 2 - u t , 

Let 

V(x,t) = \\f\\ L °°(R") +t\\g\\ L °°(Q T )- 

Then 

L[u {m) ~v]= g(x, t) + | \g\\ L <~(Q T ) > 0. 
On the parabolic boundary of B R (0) x [0,T], 

u (m) _ v < Q 

By the comparison principle, we can see that 

u (m) (x,t) < v(x,t) < M. 

The lower bound of can be derived similarly. In general, we use the standard 
approximation (see page 715 of [5], for example) to derive the desired estimate. 
To prove the second conclusion, we note that 
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R n 



u {m \x,t)dx- / f{x)dx 



R" JO 



g(x, r)drdx. 



It follows that 



ti 



(m) I 



Li(R") ^ II/IUh-r 2 ) + \\9\\l\q t )- 



For any fixed smooth function 4>(x), we have 

f A{u^) m <i){x)dx = f {u {m) ) m A<P{x)dx 

JR n JR n 

< CM m " 1 ||« (m) || i i( fl n) — >• asm^oo. 



It follows that 

/ u (oo \x,t)(P(x)dx 

JR n 

Thus, 



R n 



f(x)(j)(x)dx 



R n 



g(x,r)dr 



(f)(x)dx. 



u {qo) {x 



,t) = f( x ) + / g(x,r)di 
Jo 



for a.e. x G R r \t > 0. 

Q.E.D. 

From Lemma 4.1, we know that 

u {rn) ig 

uniformly bounded. Then for any fixed 
t > every sequence u^ m '\x, t) with m' — > oo has a subsequece i^" 1 ") with to" — > oo 
such that 

u (m") _^ u (oo)^ a . j ^ W e a kly in * - L°°(Q). 

In general i) may not be unique. However, we show that there is a unique 

limit for u^°°\x, t) as t — > 0+. 

Lemma 4.2: Under the assumption H(4.1) the following estimate holds: 



\u t \\ LHRn) < 1 \\g\\ L i {R n ) + 



Proof: Let u(x,t) = A m - 1 n(x, At). Then it is easy to see that u solves the following 
problem: 

ut - Au m = X^g{x,t), (x, t) G R n x (0, oo), 
0) = X^f(x), x e R n 

By ^-stability for the porous medium equation ([4]) we know that 



-i 



\L l {R n )\, 



where C depends only on n, but not on m. By taking limit as A — > 1+, we obtain 

1 



m 



< 



m 



C 



m 



j\\g\\V(R") + j-H/llLiOR"), 



which yields the desired estimate. 
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Lemma 4.3: The limit solution satisfies 

0<u (oo) <l mQ T . 

Proof: For any 5 > 0, we multiply the equation (4.1) by u = [u — (1 + b~)} + to obtain 
that 

/ \u\ 2 dx + [ T I miT-^Vu^dxdt 

JR n JO JR n 



<[ \u(x,0)\ 2 dx+ [ T [ \g\ 2 dxdt<C, 

JR n JO JR n 



]2 dxdt < —. 



where C is independent of m. 
It follows that 

r i iv«i s 

Jo JR n m 
Fatou's lemma implies 

\{(x,t) > 1 + S}\ = 0, 

which yields the desired result since 5 is arbitrary. 

Q.E.D. 

From Lemma 4.3, we see that there must exist an initial collapsing instantly if 
||/||l°°(_R") > 1- The external force should play no role in this collapsing process. 
The basic strategy is to separate the initial collapsing and the evolution processes. 
H(4.2): Let L =||/| | L =o (fln) >l. 

Theorem 4.4: Under the assumptions H(4-l)-(4-%), 

t) has a unique limit as 

m — > oo and t — > 0+. Moreover, 

u (co \x,0+) = v {oo \x), xeR n , 

where (x) is the limit function of the solution to (4-l)-(4-2) with g(x,t) = 0. 
Proof: Let be the solution of (4.1)-(4.2) with g(x,t) = 0. Let <f>(x) be any 
bounded function with compact support. Let (m f ,t m i) and (m",t m » be any two 
sequences pairs which converge to (oo, 0+), respectively, as m — > oo. For any bounded 
function 

| / <f>\u^'\x,t ml )-v^"\x,t m// )]dx\ 

jR n 1 J 

f f^m 1 ( t\ /*^m" ( ff\ 

l°°(r™) [ \ut n, \dr + / \vl m 'dr]dx 

jR n Jo Jo 

max{t m ,,t m „} r C C 



< C\\(p\\ L co( Rn) 

->o, 



rmax\t m , ,t m „ i r (J 

/ / \g\dxdT-\ - 

Jo JR n m' — 1 



as m',m" — > oo. Let the limit function be denoted by 

From [5], we know that v^ m \x,t) converges to a unique limit v^°°\x) as m — > oo. 
Particularly, we see that 

v {rn "\x,t m ») -> v (oo) (^0+) = v (oo) (x) 
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as m — > oo. 
It follows that 



(f)u < - oo) (x,0+)dx = / (pv {QO \x,0+)dx. 



X 



a.e.x e R n . 



Since <p(x) is arbitrary, we see that 

u (co) (x, 0+) = v {co) (x, 0+) = v (co) ( 

Q.E.D. 

Remark 4.1: Physically, Theorem 4.3 means that the collapsing process happens 
instantly before the force term g(x, t) affects the evolution process. 

To see the profile of the limit solution for t > 0, we may start with u^°°\x, 0+) as 
an initial value. Therefore, without loss of generality, from now on we assume that 

L = ||/||l°°(.r») < l- 



Lemma 4.5 Assume that 1 1/| \l°°(r™) < 1- Then 



m — 1 v ' 



m—l 



is uniformly bounded. 

Proof: With loss of generality, we may assume that V m is smooth and suppV m C 
.Br(O) for all t E [0,T]. A direct calculation shows that V m satisfies 

V m , t = (m- l)V m AV m + \VV m \ 2 + V£g(x, t),xeB R ,te (0, T], 
V m (x,t) = 0, (x,t) edB R (0) x [0,T], 



V m (x,0) 



m 



m—l 



f(x) m ~\x e K 



where (3 = < 1- If K« attains positive maximum M > 1, then the coefficient 
of V m , ?_0 , has an upper bound which is independent of m. The strong maximum 
principle implies that V m can not take a positive maximum which is greater than 1 
in B R (0) x (0,T]. It follows that 



m 



V m < maxfl, max V m (x, 0)} < 

b r (o) m-l 



< 2. 



Q.E.D. 

Lemma 4. 6: M < - m - ) has a uniform compact support. 

Proof: By Lemma 4.5, V m is uniformly bounded. Now consider the special Baren- 
blatt solution 



e - 



U — t "(™-l)+2 

where (s)+ = max{s, 0} and 

M (2(n{m- 1) + 2 



m — l 



2[n(m - 1) + 2] 



X t n(m-l) + 2 



1 

m — l 



m — l 



n(m-l) + 2 
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Assume that suppf, suppg C B R (0) for some large R and for every t G [0, T\. 
By choosing M sufficiently large and using the comparison principle in the domain 
(R n \B R (0)) x (0,T], we obtain the desired result. 

Q.E.D. 

H(4.3): Assume that f(x) G C 1 (i? n ) is radially symmetric with compact support 
and f r < in i? n \{0} f] suppf. Assume that g(x,t) G C 2 ' l (R n x [0,T]) is radially 
symmetric in x and g r < in i?™\{0}) f| suppg for each t G [0, T]. Moreover, 

Af m + g(x,0)>0. 

Lemma 4.7: Under the assumptions H(4.1)-(4.3), u^> < and uf 11 ^ > 0. 
Proof: This can be derived by a simple application of the maximum principle. 

Theorem 4.8: Under the assumptions H(4.1) and H(4.3), the limit solution is 
uniquely determined and 

1, in {(x,t) : w(x,t) > 0}, 
f(x) + Jo g(x,r)dr, in {(x,t) : w(x,t) = 0}, 

where w(x,t) is the unique solution of the variational inequality 

-Aw>f(x)+ f g(x,r)dr-l, (4.3) 
Jo 

w(x,t)>0, (4.4) 

w(x, t) [Aw + f(x) + C g(x, r)dr - 1] = (4.5) 
Jo 

for a.e. (x,t) G R n x (0,T]. 

Proof: Since the proof closely follows the argument of [5] (Theorem 5.2-5.3), we only 
give a brief outline here. First of all we recall that 1 1/| \l°°(r™) < 1- Since f r < 0, 
we may further assume that f(x) = 1 on a small ball, say, Bg(0) for some small 
8 > 0. Since g(x, t) > 0, by the comparison principle (see [5] Theorem 3.2) we see 
that ii(°°)(a;,*) = 1 for a.e. x G B s and t G [0,T]. 
Step 1: For each t G [0,T], define 

N t = {x:f(x)+ f g(x,r)dr=l}. 
Jo 

By the comparison principle, it is clear that N tl C N t2 if ti <t 2 . We claim that 

u (°°\x,t) = l,for a.e. all x G N t . 

Indeed, we replace / and g by f £ = min{\ — e, /} and g — e, respectively, for small 
e > 0. Let i4 m ) and -u^ 00 ^ 1 be the solutions of (4.1)-(4.2) corresponding to f £ and g £ , 
respectively. Then Lemma 4.1 yields 

4°°\x,t) = f £ (x)+ f g £ (x,T)dr 

J 



u {oo \x,t) = 
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for a.e. all x G R n . On the other hand, the comparison principle implies 

u {m) {x,t) > 4 rn \x,t). 

It follows that 

u {QO \x,t) > f £ {x) + f\g-e)dr. 

Jo 

It follows by the definition of N t that 

u (°°\x,t) = 1, ase^O, 

for a.e. x G N t . 
Step 2: Let 

S t = {x : w (oo) (:r,t) < 1}, 

S% = {x G St '■ St has Lebesgue density 1 at x}. 

Assume that for any xq G S% and to £ (0,2""] with u^°°\xo,to) < 6o < 1. Then for 
any 9 G (#o, 1), 5 > 0, there exists an mo such that 

inf W ( m ')(:r,t) < 0, 
xeQi(xo,to) 

if m' > m . 

The proof is almost the same as the proof of Theorem 3.3 in [5]. We skip it here. 
Step 3: If x G S* with u (oo) (x , to) < < 1, then 



«(°°>(a;,t) = /(a;)+ f g{x,r)dr 
Jo 



in the (i?"\ J B ro (0)) x [0,T ], where r = \x \. 

This step can be proved by using Step 3 and Lemma 4.7. Indeed, by Step 2 we 
know that for any 9 G (6> , 1) and any 5 > there exists m such that 

inf u {m '\x,t Q ) < 9 

Bs(x ) 

if m! > m . Suppose that the minimum of u^ m '^ attains at a point (x*,t ) in Bs(xq). 
Since < and u[ m ^ > 0, x* must occur on the boundary of B$(x ) and, if 
m! > m , then 

w ( m ')(x,t) < e 

for all (x, t) G (i? n \S ro+ 5(0)) x [0,t ], where r = |x |. It follows by Lemma 4.1 and 
Lemma 4.7 that 

«(°°>(M) = /(a;) + fg{x,T)dT 
Jo 

for all t) G (R n \B ro+ s(0)) x [0, £ ]- Finally, since <5 is arbitrary, we conclude the 
desired result. 
Step 4: 

u {oo \x,t) = XN t + [f(x) + ^ g(x,T)dT]x R ™\Nf 
Indeed, by Step 1 to Step 3, there exists a function < f(x, t) < 1 such that 

U (00) (x,t) = XN t + f(x,t)xnn\Nf 
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If at a point (x , t ) with 

rto 

f(x ) + / g(x,r)dr < 1, 
Jo 

then by Lemma 4.1 and Lemma 4.7, 

u {QO \x,t) = fix) + f g(x,r)dr 
Jo 

for all x G R n \B ro (0). It follows that 

!{x,t) = f(x) + / g(x,r)dr. 
Jo 

Step 5: 

iVf = {x : £) > 0}, 

where w(x,t) is the solution of the variational inequality (4.3)-(4.5). 

First of all, by Lemma 4.6 we know that has compact support. Define 

w(x ,t) = J Rn (f(x) + j\{x,r) - u (oo) (^*)) T *odx, 

where T XQ is the fundamental solution with singularity at xq. Note that w(xo,t) is 
well defined since has compact support. It is clear that 

— Aw = fix) + / gix, r)dr — 1 
Jo 

in N t (since u^°°\x, t) = 1 on N t ). 
On the other hand, 

w^ m \x ,t) = ( [u^\x,t)~ f{x)- f g{x,r)dr]Y xo dx 

jR n Jo 

= f I A(u^) m T Xo dxdT 
Jo Jn n 

= [\ U W) m (xO,T)dT. 

Jo 

Clearly, w^ m \xQ,t) > for any x G R n and t G [0,T]. By Step 4, we see that if 
x G St, then there exists a small ball Bs'{xq) C St. This implies 

w (m,) (^o,r) < 9 < 1 

for all r G [0, t], provided that m' is sufficiently large. Consequently, w^ m \xo,t) — > 
as m — > oo. Thus, 

> 0, for a.e. x eR n ,te [0,T], 
w(xo, £) = 0, if xo G iV t , 
i.e. w(x,t) is a solution of the variational inequality (4.3)-(4.5). 

Q.E.D. 

5. The Profile of the Current Density for Limit Solution 
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Throughout this section we assume that 

U = {h 1 ,h 2 ,0}. 

We shall use the results from the previous section to derive the profile of the current 
density. 

H(5.1): Assume that 

L = ||V x H || L -(ij2) > 1. 

H(5.2): Assume that F and H satisfy the assumptions H(4.1)-H(4.2), where u 0i and 
Qi are defined as in section 2. 
First of all, since 

VH = /ioixi + ^02x2 — 0) 

we see by the definition that 

u Q1 (x) = u 02 (x) = h 02x2 (x) - h 01xi (x). 
Similarly, since V x F = one can find 

gi(x,t) = 92(x,t), inQ T . 
Define the current densities J p and J as follows: 

J p (x,t) = |V x H (p) |, J(x,t) = |V x H (oo) |. 
Note that J is well-defined since V x H^ 00 - 1 is well-defined. 

Theorem 5.1: Under the assumptions H(5.1)-(5.2) the current density J(x, 0+) = 
1 if x e So and J(x,0+) = |V x H | if x e N , where N = R 2 \S , S is the 
noncoincident set of a variational inequality 

-Aw >g-l,w>0, (Aw + g t - l)w = 0. 
Proof: From Theorem 4.3, we know that 

u ( f°\x,0+) =v^\x). 
On the other hand, we know that 

u (p) = au (p) = | V x h (p)|. 

It follows 

J p (x, i) — > J(x) 

as p — > +oo. Thus, in the supercounductor region, denoted by So, the current density 

J(x,0+) = ^ (oo) (x) = |V x H |, a; G i? 2 
In the normal conductor region, denoted by N$, the current density 

J(x,0+) = v {oo) (x) = I. 
Moreover, S is the noncoincident set of a variational inequality (see [5]): 
—Aw > ua — 1, w > 0, (Aw + u 0i — l)w = 0. 

Q.E.D. 
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Next we assume that 

||V x H ||L°°(i?2) < I. 

H(5.2)': Assume that H and F satisfy the assumptions H(4.3), where u 0i and g 0i 
are defined the same as in Section 2. 

Theorem 5.2: Under the conditions H(5.1)-H(5.2) ) the current density J(x,t) = 
|V x H(°°)| has the following structure: 

f i,ifxeN t ,t>0, 

1 ' } \ f{x) + fo9(x,r)dr, if x G S t ,t > 0. 
where N t = R 2 \S t and S t is the noncoincident set of the variational inequality 

w(x, t) > 0, — Aw > u 0i + / g,idr — 1, 

Jo 

w(Aw + u i + / gidr - 1) = 0. 
Jo 

Proof: From Theorem 4.8, we see that 

J(x,t)=ut\x,t). 
The rest of the proof follows the same procedure as that of Theorem 5.1. 

Q.E.D. 

Next we show the limit solution of (1.1)-(1.3) satisfies a degenerate evolution sys- 
tem. 

Theorem 5.3: There exists a nonnegative, bounded and measurable function a(x,t) 
such that the limit solution H^°°^ satisfies 

H^ oo) - V x [a(x, t)V x H(°°>] = F, (x, t) G Q T , 

in the sense of distribution. Moreover, 

suppa C N t for a.e. every t G (0, T], 

Proof: The proof is similar to [9]. Define 

A p = |V x H| P ~ 2 VH. 
Then, for any V G H{0, T; B d ) with suppV G B R (0) for each t G (0, T], we have 

[ T [ [-H (p) • V t + A p • VV]dxdt 
Jo Jb r (o) 

Jo Jb 

By Lemma 4.6, we know that there exists a subsequence such that 
A Pfc -> A, weakly-* in L°°(R 2 x (0, T)). 



rT 

F • Vdxdt. 

IB R (0) 
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It follows that 



[ T [ [-H (oo) • V t + A • V x V]dxdt 
Jo Jb r (o) 



F • Vdxdt. 

o Jb r (o) 



On the other hand, we take V = H^ p ^ as a test vector to obtain 

1 [ \H^\ 2 dx+ [ T [ \VK {p) \ p dxdt 

2 Jb r (0) JO JB R (0) 



'Bfl(O) Jo JB R (0) 

rT r , . rT 



f f F • (H&tyfodf -> I I F • (H^)dxdt. 
Jo Jb r (o) Jo Jb r (o) 



>B R (0) Jo JB R (0) 

It follows that 

r'l 

Ib r (o) 1 1 Jo Jb r (o) 



( T ( |V x U. {p) \ p dxdt -> ( T { A • (V x B_)dxdt 

Jo Jb r (o) Jo Jb r (o) 



as p — > oo. 

By Fatou's lemma, we have 

i 7^(0)' ' ~ JPk Jo JB R (0) Pkl 

< lirn^ {j* J b (q) |V x U^\ Pk \dxd?j Pk [T\B R (0)\% 

= !! A • (V x H^dxdt. 

Jo Jb r (o) 

It follows that 

|A| = A - (V x H (oo) ), 

since |V x H^°°^| < I. Consequently, there exists a nonnegative, bounded function 
a(x, t) such that 

A(x,t) = a(x,t)V x H (oo) . 
Moreover, if V x exists and |V x H<°°)| < 1, then a(x,t) = 0. 

Q.E.D. 
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